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ABSTRACT 

We use cosmological simulations in order to study the effects of supernova (SN) feed- 
back on the formation of a Milky Way-type galaxy of virial mass ^ 10^^ /I'^M©. 
CN I We analyse a set of simulations run with the code described by Scannapieco et al. 

(2005, 2006), where we have tested our star formation and feedback prescription using 

' (-H '' isolated galaxy models. Here we extend this work by simulating the formation of a 

O „ galaxy in its proper cosmological framework, focusing on the ability of the model to 

I ■ form a disk-like structure in rotational support. We find that SN feedback plays a fun- 

O ' damental role in the evolution of the simulated galaxy, efficiently regulating the star 

formation activity, pressurizing the gas and generating mass-loaded galactic winds. 
^ ' These processes affect several galactic properties such as final stellar mass, morphol- 

ogy, angular momentum, chemical properties, and final gas and baryon fractions. In 
particular, we find that our model is able to reproduce extended disk components with 
Cn ' high specific angular momentum and a significant fraction of young stars. The galaxies 

^ I are also found to have significant spheroids composed almost entirely of stars formed 

at early times. We find that most combinations of the input parameters yield disk-like 
components, although with different sizes and thicknesses, indicating that the code 
can form disks without fine-tuning the implemented physics. We also show how our 



in 



model scales to smaller systems. By analysing simulations of virial masses 10^ h ^Mq 



^T ' and 10^" h ^M©, we find that the smaller the galaxy, the stronger the SN feedback 

O ; effects. 

OO 

f^ . Key words: galaxies: formation - evolution - abundances - cosmology: theory - meth- 

ods: numerical 
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1 INTRODUCTION lar disks are thickened or destroyed by rapid fluctuations in 

the gravitational potential such as those produced by the 
The formation of spiral galaxies is currently a subject of ^- r ^ n-^ i • ^rn ^i c ^ _l -i -i^,ir, r^ ■ 

, , , accretion oi satellite galaxies ( loth & Ostriker 1992; (jumn, 

considerable interest m cosmology. From a theoretical point ^t • . o i-i h -i^^o ^r , o tm •. ^nnn\ mi 

, . , ^ . r ,. , n • . „ , Hernquist & FuUagar 1993; Velazquez & White 1999). The 

oi view, the tormation ot disk coniigurations is well under- r .i . i. i . . i . i 

presence or tnm disks m many present-day spirals suggests 
stood as a result ot specinc angular momentum conservation , , , , , , , i i- i • ^ j_- ^ ■ ^ 

, , p that they had a relatively quiet mass accretion history, at 

during the collapse and relaxation oi gas withm dark matter , . . . i 

i63,st in tiiG rscsntj "Deist 
potential wells (Fall & Efstathiou 1980; Mo, Mao & White 

1998). Dark matter and baryons aquire similar speciflc angu- ^-r .,.,,. , , . , . , i 

' . , , , r ,, /TT , JNumerical simulations have become an important tool 

lar momenta irom tidal torques before collapse (Hoyle 1953; r ^ ^ ■ i r ^- i ^ ^i • i -i-^ ^ j -i. 

TD ui men ijri,-4- iaoA\-D -J J ^-u ^ ^ '' .' for studying galaxy formation due to their ability to describe 

Peebles 1969; White 1984). Provided the gas conserves this ,, ,. , , , ,. r j i ,, j i_ 

' ' ° the nonlinear, coupled evolution ot dark matter and baryons 

angular momentum, its collapse leads to a disk-like configu- ... i • i j_ ^ tt ii,- i ^ j ^ i 

° J ±' ° m its cosmological context. However, this has turned out to 

ration at the centre of its dark matter halo. If stars form from , , i i n • ii i r J^^ i^ 

be an extremely challenging problem because ot the need to 
this rotationally supported gas, then spiral disks like those j -i, ■ i^ i ii. i i j • r i 

■' '^^^ ° ' '^ describe simultaneously the large-scale dynamics ot galaxy 

observed can develop. However, centrifugally supported stel- , , , , , n i i ^ j ^ i c 

'^ ' a J ff assembly and the small-scale processes related to star torma- 

tion and evolution. State-of-the-art simulations have imple- 
mented schematic recipes for radiative cooling, star forma- 
* E-mail: ceciliaOmpa-garching. mpg.de tion, feedback from supernovae (SN), stellar winds, black 
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hole formation and active galactic nuclei. Although these 
are successful in many respects, they have had difficulty re- 
producing extended disk-dominated galaxies in typical dark 
matter haloes. The main problem is so-called catastrophic 
angular momentum loss, which arises when the baryonic 
components transfer most of their angular momentum to 
the dark matter during condensation into the central re- 
gions, forming dominant spheroidal stellar systems and disks 
that are too small compared to observations (e.g. Navarro 
& Benz 1991; Navarro & White 1994). 

The formation of disks is heavily influenced by details 
of the gas cooling process. Disks that are less centrally con- 
centrated and with higher specific angular momentum are 
obtained if cooling is artificially suppressed during the as- 
sembly of dark matter haloes (Weil, Eke & Efstathiou 1998; 
Eke, Efstathiou & Wright 2000). Many recent studies have 
attempted to prevent the early cooling of gas and the loss 
of angular momentum, some of them successfully reproduc- 
ing disks as big as observed spirals (e.g. Abadi et al. 2003; 
Governato et al. 2004; Robertson et al. 2004; Okamoto et 
al. 2005; Governato et al. 2007). These studies have pointed 
out the importance of stellar feedback as a key process pro- 
viding pressure support to the gas and preventing its rapid 
collapse within the haloes, thereby regulating star formation 
activity. However, it seems that fine-tuning of the models is 
required in order to reproduce the observed properties of 
disk galaxies in a hierarchical scenario, and a detailed un- 
derstanding of disk formation remains elusive. Insufficient 
numerical resolution or a poor description of the interstellar 
medium (which influences the efficiency of the stellar feed- 
back) have also been proposed as causes for the inability 
to reproduce convincing disk galaxies in simulations (e.g. 
Okamoto et al. 2003; Governato et al. 2004). 

Over the years, many authors have included star forma- 
tion and feedback recipes both in mesh-based and in SPH 
codes in order to study the effects of SN feedback on galaxy 
properties such as star formation histories, metallicity evolu- 
tion, generation of galactic outflows, abundance of satellites 
and galaxy scaling relations (e.g. Katz & Gunn 1991; Gen 
& Ostriker 1992; Navarro & White 1993; Metzler & Evrard 
1994; Steinmetz & Miiller 1995; Gerritsen & Icke 1997; Yepes 
et al. 1997; Gen & Ostriker 1999; Sommer-Larsen, Gelato & 
Vedel 1999; Thacker & Gouchman 2000; Kay et al. 2002; 
Semelin & Gombes 2002; Abadi et al. 2003; Marri & White 
2003; Sommer-Larsen, Gotz & Portinari 2003; Springel & 
Hernquist 2003; Brook et al. 2004; Okamoto et al. 2005; Op- 
penheimer & Dave 2006; Stinson et al. 2006; Governato et al. 
2007; Dalla Vecchia & Schaye 2008; Dubois & Teyssier 2008; 
Finlator & Dave 2008) . All these feedback implementations 
require arbitrary and ad hoc elements because the relevant 
dynamical processes cannot be resolved by the simulations 
and so have to be introduced as "sub-grid" models. As a 
result, their validity is uncertain and the conclusions to be 
drawn from their successes and failures are correspondingly 
unclear. A further difficulty is introduced by numerical lim- 
itations which have so far prevented the simulation at high 
resolution of a large and representative sample of galaxies 
in their proper cosmological setting; typically, the evolution 
of a small number of selected objects is studied in detail, 
often chosen as those most likely to produce disk galaxies. 
It is thus uncertain whether existing models can successfully 
reproduce the variety of disk galaxies observed in the Uni- 



verse. In particular, no AGDM simulation has yet managed 
to produce a late-type spiral like M33 or MlOl. 

In Scannapieco et al. (2005, 2006), we have presented a 
new model for the effects of SN feedback in cosmological sim- 
ulations. This model was implemented within the Tree-PM 
SPH code GADGET-2 (Springel & Hernquist 2002; Springel 
2005), and is described in detail in Scannapieco et al. (2005, 
2006). Our model includes a scheme to treat a multiphase 
interstellar medium which allows the representation of a co- 
spatial mixture of cold and hot interstellar medium (ISM) 
components without introducing scale-dependent parame- 
ters. Energy and chemical feedback by Type II and Type la 
SNe, as well as metal-dependent gas cooling (Sutherland & 
Dopita 1993), are also considered. The model is intended to 
be as simple as possible and to introduce the fewest input 
parameters possible. It was specially designed to avoid the 
inclusion of mass-scale dependent parameters. This makes 
it suitable for running cosmological simulations where, at 
any given time, systems of different masses are forming si- 
multaneously. In Scannapieco et al. (2005, 2006) we tested 
the code by studying the formation of disk galaxies of dif- 
ferent mass from idealized initial conditions. Our scheme 
can reproduce several important effects in galaxy evolution: 
substantial metal-enhanced galactic winds are produced, re- 
distributing mass and metals and generating a self-regulated 
cycle for star formation. The galactic winds are generated 
naturally in our model when the multiphase gas model and 
SN feedback are both included. This is an important result 
which allows us to consistently study the formation of galax- 
ies together with the chemical enrichment and the redistri- 
bution of matter within the galaxies. Potentially it should 
also allow investigation of other aspects such as the enrich- 
ment of the intergalactic medium. Our model reproduces 
the expected dependence on galaxy mass: while star forma- 
tion is suppressed at most by a factor of a few in massive 
galaxies, in low-mass systems the effects can be much larger, 
giving star formation an episodic, bursty character. In this 
paper, we extend the work of Scannapieco et al. (2005, 2006) 
by studying the formation of disk galaxies in their proper 
cosmological context. This allows a consistent treatment of 
processes such as mergers, interactions, inffows, tidal strip- 
ping, etc. We thus carry out higher resolution simulations 
from realistic initial conditions and follow the hierarchical 
growth of a galaxy similar in size and spin to our Milky 
Way. We focus on the effects of SN feedback on star forma- 
tion and on galaxy morphology, emphasizing the dependence 
of the predictions on input parameters. 

Our paper is organized as follows. In Section [5] we de- 
scribe the simulation code and list the characteristics of our 
simulations, brieffy describing the inputs which are relevant 
for this work. In Section [3] we explore the effects of SN feed- 
back on star formation fsubsection l3.1|) . galaxy morphology 
(subection 13. 2[) and evolution of specific angular momen- 
tum (subsection 13. 3p . We also study the formation process 
of the disk components found in our simulated galaxies (sub- 
section 13. 4p and the resulting chemical properties (subsec- 
tion [331 ■ In Section m we investigate how SN feedback pro- 
ceeds for smaller galaxies. Finally, we give our conclusions 
in Section [5l 
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2 THE SIMULATION CODE 

The simulations in this paper were run with the code de- 
scribed in Scannapieco et al. (2005, 2006), which is an exten- 
sion of the parallel Tree-PM SPH code GADGET-2 (Springel 
& Hernquist 2002; Springel 2005). The code includes star 
formation, metal-dependent cooling, chemical enrichment 
and energy feedback by Type II and Type la Supernovae, 
in addition to a multiphase treatment for the gas compo- 
nent. At the end of this Section, we describe those input pa- 
rameters of our model that are particularly relevant for the 
present study, and we refer the interested reader to Scanna- 
pieco et al. (2005) for details on the implementation of star 
formation, chemical enrichment and metal-dependent cool- 
ing, and to Scannapieco et al. (2006) for the description of 
the multiphase gas and energy feedback models. Note that 
our multiphase treatment and our star formation and feed- 
back models differ substantially from those of Springel & 
Hernquist (2003) although we do adopt their treatment of 
the UV background. 

We focus on a study of the formation of disk galaxies 
in their cosmological context. For this purpose we simulate 
a system with z — halo mass ~ 10^^ h~^ Mq, extracted 
from a large cosmological simulation and resimulated with 
improved resolution. Our simulation runs from z = 38 to 
z = and adopts a ACDM Universe with the following 
cosmological parameters: Qa = 0.7, Qm ~ 0.3, flh ~ 0.04, a 
normalization of the power spectrum of erg = 0.9 and Ho = 
100 h km s~^ Mpc~^ with h — 0.7. The particle mass is 
1.6 X 10^ for dark matter and 2.4 x 10® h~^ M© for baryonic 
particles, and we use a maximum gravitational softening of 
0.8 h~^ kpc for gas, dark matter and star particles. 

At z — the halo of our galaxy has spin parameter 
of A = 0.03 and it contains ~ 1.2 x 10^ dark matter and 
~ 1.5 X 10^ baryonic particles within the virial radius. It 
was selected to have no major mergers since z = 1 in order 
to give time for a disk to form. Because the selected halo 
is similar in mass and spin to the Milky Way and has a 
relatively quiet merger history in the recent past, we thus 
expect this halo to host a disk galaxy. However, given our 
incomplete knowledge of galaxy formation processes, this 
expectation is not necessarily fully correct. 

We have run a series of simulations in order to study 
the dependence of the results on assumed input parame- 
ters and to see if the code can produce disk systems. Be- 
cause of this focus, we decided to vary only the parame- 
ters related to star formation and energy feedback, and we 
fixed those related to chemical enrichment and phase de- 
coupling. In all our simulations we have assumed a SNIa 
rate of 0.3 relative to SNII and a SNIa lifetime uniformally 
distributed in the range [0.1, 1] Gyr. We have adopted the in- 
stantaneous recycling approximation for SNII, we have used 
metal-dependent cooling (Sutherland & Dopita 1993) and 
we have adopted a Salpeter Initial Mass Function. Chemical 
yields from Woosley & Weaver (1995) and Thielemann et 
al. (1993) have been used for SNII and SNIa, respectively. 
The multiphase gas model of Scannapieco et al. (2006) has 
been turned on for all experiments. This model decouples 
particles with dissimilar properties, preventing them from 
being SPH neighbours. In all the simulations presented in 
this paper, we have used a = 50 for the decoupling parame- 



ter which sets the difference in entropy required for particles 
to be decoupled (see Scannapieco et al. 2006 for details). 

For the star formation and feedback models, there are 
three relevant input parameters: a) the star formation ef- 
ficiency c, b) the energy _Bsn released per SlSo and c) the 
so-called feedback parameter ec- In our model, we assume 
that gas particles are eligible for star formation if they are 
denser than a critical value (p > pcrit = 7 x 10~^®g cm~'^ 
where p denotes gas density) and lie in a convergent flow. 
For these particles, we assume a star formation rate (SFR) 
per unit volume equal to 



p* =c- 



Tdyn 



(1) 



where c is the star formation efficiency, p« denotes stellar 
density, and Tdyn = l/V4vrGp is the local dynamical time of 
the particle. The star formation efficiency then sets the effi- 
ciency at which stars are formed or, equivalently, the typical 
time-scale of the star formation process. Once star particles 
are formed, they are assumed to give rise to both SNII and 
SNIa explosions (see Scannapieco et al. 2005 for details on 
the implementation of star formation and of SN explosions) . 
Each SN event will release an amount of energy, given by the 
number of exploding stars within the stellar particle multi- 
plied by the energy per SN, -Esn. Finally, a feedback pa- 
rameter is introduced in our model when implementing the 
transfer of SN energy to the gas. This is based on a separa- 
tion of gas particles into two phases called hot and cold (split 
by assumed thresholds in density and temperature). When 
SN explosions take place, the energy is distributed as ther- 
mal feedback in a fixed proportion to these two phases, the 
cold phase receives a fraction e^ and the remaining fraction 
1 — Ec is distributed to the hot neighbours. While hot par- 
ticles receive the energy instantaneously, each cold particle 
instead accumulates the SN energy until it is large enough 
to raise its entropy to the level of its local hot neighbours. At 
each time we define for each cold particle a local hot phase 
and compare the average entropy of the identified hot neigh- 
bours with the entropy that the cold particle would have if 
it absorbed the accumulated energy. This gives the code the 
flexibility to accomodate itself to local conditions, since dif- 
ferent systems will have different typical hot-phase entropies 
and these will change with time. The parameter ec is there- 
fore related to the mode of feedback in a non-trivial way 
(Scannapieco et al. 2006). Together with the values adopted 
for the star formation efficiency and energy released per SN, 
the feedback parameter determines the details of the evolu- 
tion of the galaxy. 

A list of the simulations we have carried out is shown 
in Table [l] They are divided into three series: F, E and 
C refering to the parameter which is varied in each series. 
Simulations F were run with different feedback parameter 
(ec), simulations E with different energy per SN (-Esn) and 
simulations C with different star formation efficiency (c). 
Note that F-0.5 is a standard case which belongs to all three 
groups of simulations. We have also run a simulation without 
energy feedback (although including chemical enrichment. 



^ In Scannapieco et al. (2006) we used a parameter e^ which 
denotes the fraction of the SN energy assumed to be lost by ra- 
diation. In practice, this is equivalent to assume a lower value for 
Esn- In this work, we decided to use Esn instead of e^. 
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Table 1. Principal characteristics of our simulations: fraction of energy (and metals) distributed into the cold phase (ec), energy per SN 
(EgN in units of 10^^ ergs) and star formation efficiency (c). We also list dark matter, gas and stellar mass in units of 10^'^ /i^^Mq, gas 
fraction and baryonic fraction both within the virial radius (superscripts '200') and within twice the optical radius (superscripts '2ropt'). 
We show the optical radius for the different simulations, in units of h^^kpc. The fraction of stellar mass formed at z < 1 is also shown. 
All quantities are computed at z = 0. 
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metal-dependent cooling and the multiphase gas treatment) 
in order to highlight the effects of this process. This simula- 
tion is referred to as NF. 



3 IMPACT OF SN FEEDBACK ON THE 
FORMATION OF GALAXIES 

In this Section we use our simulations to investigate how our 
model of SN feedback affects the formation of a galaxy and, 
in particular, if this process affects the formation of a disk- 
like structure. We also analyse the dependence of our results 
on the values of the relevant input parameters. We will show 
that SN feedback plays a fundamental role, regulating the 
star formation rate, pressurizing the gas and allowing its 
late collapse to a disk structure. In the next Section, we 
investigate how SN feedback affects the evolution of two 
galaxies with two and three orders of magnitude smaller 
mass. 

In Table [T] we list a variety of quantities for our sim- 
ulations: the masses of dark matter, gas and stars (A^dm, 
Mgas and Afstar), the gas fraction /g (i.e. the ratio between 
the gaseous and baryonic mass), and the baryon fraction 
/bar (i-e- ratio between baryonic and total mass). All these 
quantities are computed at z = and refer either to the 
region within r2oo, the radius enclosing a mean density 200 
times the critical value, or to the central regions (r < 2ropt, 
where ropt is defined as the radius which encloses 83 per- 
cent of the baryonic mass of the central object). Note that 
our haloes are formed by a central system (refered to as the 
central galaxy) and a series of satellites. Hence, quantities 
inside the virial radius are affected by the properties of the 
satellites, whereas those inside twice the optical radius give 
a cleaner estimation of the characteristics of the central sys- 
tem. In the last column we show /^^^, which denotes the 
total stellar mass formed at z < 1 divided by the final stellar 
mass at z = (calculated within r < 2ropt). 



3.1 Effects on the star formation process 

In Fig. [l]we show star formation rates (SFR) for our sim- 
ulations as a function of time. We have split our plots to 



highlight the dependence of the results on each of the in- 
put parameters. The top panel shows the SFRs for different 
Ec values (F-0.3, F-0.5 and F-0.9) and for the run without 
energy feedback (NF). SN feedback produces a substantial 
decrease in the SFR, especially at early times. In the no- 
feedback case, the gas collapses and concentrates at the cen- 
tre of the potential well, producing a strong starburst which 
only fades when the amount of available gas decreases. Since 
there is no mechanism to expel gas in this case, the SFR be- 
haves smoothly following the gas collapse. On the contrary, 
when SN feedback is considered, the level of SFR is signif- 
icantly lower at early times and its bursty behaviour indi- 
cates the self-regulation of star formation activity produced 
by the injection of SN energy into the interstellar medium. 
In this case, each starburst is followed by a depression in the 
SFR because the expansion of surrounding gas decreases the 
amount of material available for star formation. 

Our results are sensitive to the adopted value for the 
feedback parameter tc. For larger e^ values, the cold gas 
receives more energy, and hence we might expect stronger 
regulation effects. Note that energy injection to both hot and 
cold phases affects the evolution of the system. However, the 
cold gas is the fuel for star formation, and the effects of SN 
feedback on the SFR are more sensitive to energy deposition 
in the cold phase. Our results indicate that increasing ec 
reduces the SFR at early times but can result in more gas 
being available for star formation at later phases. 

From Table [T] we see that the no-feedback simulation 
(NF) forms 24.6 x 10^" /I'^M© of stars inside the virial ra- 
dius, and 14.9 x 10^° h~'^Mrs, within 2ropt. The difference 
between these two values reflects the presence of massive 
satellites which contribute signiflcantly to the stellar mass 
inside r20CLl- As a result of efficient cooling and star forma- 



^ Note that our simulated galaxies are formed by a central system 
and a series of satellites. In order to have a better estimation of 
the contribution of satellites to the stellar mass of the simulated 
galaxy, we have done a proper decomposition of the mass in the 
halo into a series of subhaloes, and calculated the stellar mass 
contributed by each of these subhaloes to the total mass. In this 
way, we find that 25 percent of the total stellar mass is contributed 
by satellites in the final state of the simulation. 
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Figure 1. Star formation rates for tlie tliree different series of 
experiments. Upper panel; results for varying feedback parameter 
(F-0.3, F-0.5 and F-0.9), as well as for the no-feedback simulation 
(NF). Middle panel: results for varying energy per SN (E-0.3, E- 
0.7, F-0.5 (-©SN = 1 X 10^1 ergs) and E-3). Lower panel; results 
for varying star formation efficiency (C-0.01, F-0.5 (c = 0.1) and 
C-0.5). Note that F-0.5 is contained in all three series. These star 
formation rates were calculated using all the stars which ended up 
within 2 Topt of the centre of the final galaxy. They thus exclude 
stars which end up in satellite galaxies. 



tion, this simulation ends up with only 22 percent of the 
baryons in the form of gas within the virial radius, and only 
5 percent in the inner regions. The inclusion of energy feed- 
back significantly changes these numbers, producing a sub- 
stantial decrease in the number of stars formed in the galaxy 
(see Table [T|. In our feedback cases, the systems are able to 
keep ~ 50 percent of the baryons within the virial radius in 
the form of gas at z = 0, and from 5 to 10 percent in the 
inner regions. Note that in the feedback cases, ~ 80 percent 
of the stellar mass within the virial radius is in the main 
galaxy, so the satellites contribute much less to the stellar 
mass at z = 0. 

The baryon fraction of the galaxies (listed in Table [1} 
is also affected by the inclusion of SN energy feedback. In 
the no-feedback case (NF), /^°° = 0.14 and fl^"' = 0.31, 
while the feedback cases have /^2r = [0.05 — 0.08] and 
/bar'"' ~ [0.14 — 0.17] which are significantly lower. This 
is a consequence of galactic winds driven by gas which is 
heated by the SN energy. 

The middle panel of Fig. [l] shows SFRs for the ex- 
periments with varying energy per SN. We compare F-0.5 
{EsN = 1), E-0.3, E-0.7 and E-3. In this case, higher _Bsn 
values translate into lower SFRs and lower final stellar and 
gas masses (see Table [!}. In these cases, the gas fraction 
within the virial radius increases with the adopted energy 
per SN, while the baryon fraction decreases. These results 
are a direct consequence of the higher energy input which 
prevents further star formation and produces stronger galac- 
tic winds. The effects of SN feedback are very clear in the 
extreme case of very high SN energy (3 x 10^^ ergs per SN) , 
where SF occurs in a series of starbursts and its level is very 
low. We will show in the next subsection that this system is 
not able to form a realistic disk. In this case, galactic winds 
are very strong, resulting in a low final stellar mass, and in 
a baryon fraction which is less than 0.05, both within the 
virial radius and in the innermost region. 

Finally, in the lower panel of Fig. [T] we show SFRs for 
the series with different star formation efficiencies (c). In the 
absence of SN feedback, a higher c would simply translate 
into a higher SFR. However, this is not the case when we 
consider the injection of SN energy into the ISM. Compar- 
ing F-0.5 (with c = 0.1) with C-0.01 and C0.5, we find that 
higher c values lead to lower SFRs and hence lower stellar 
masses at 2 = 0. This is because for a given gas distribution, 
larger c values produce more stars and so more SN feedback. 
The surrounding ISM is thus more strongly heated and re- 
acts by expanding more rapidly. This results in the effective 
suppression of further cooling and cuts off the gas supply for 
further star formation. In contrast, a lower c value results in 
a lower star formation rate and the reduced feedback heat- 
ing is unable to prevent further cooling of gas into the star- 
forming phase, with the result that over somewhat longer 
timescales more stars are formed. It is this competition be- 
tween SN heating and cooling of gas to provide fuel for star 
formation which produces the trend between c and the level 
of SFR in our simulations. For example, run C-0.01 forms 
twice as many stars as run C-0.5. 

We have shown that changing the input parameters 
of our feedback model produces changes in the SFRs, the 
stellar masses, and the gas and baryon fractions of the fi- 
nal galaxies. Despite these differences, most of our simula- 
tions reproduce approximately the Kennicutt law (Kenni- 
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cutt 1998) at 2: = 0, with the exception of our no-feedback 
model and the most extreme and least realistic feedback 
cases E-3 and C-0.5, which also fail to form disks (see next 
subsection) . 

Self regulation of star formation of the kind found in 
our models is crucial for the evolution of the galaxies. Gas 
left over from early times as a result of SN energy input 
can be maintained in a pressure supported halo which later 
collapses to produce significant star formation activity in 
the last 8 Gyr (z < 1). In the following we will see that 
this young stellar population is a major contributor to the 
formation of disks. 



3.2 Effects on galaxy morphology 

In order to see how SN feedback affects the formation of 
disks, we now examine the stellar mass distributions of our 
simulated galaxies. Fig.[2]shows edge-on stellar surface mass 
density maps at z = 0. To make them we compute the to- 
tal angular momentum of the stars in the central regions 
(r < 2 Topt), this gives a preferred direction which we define 
as the 2-axis (the angular momentum is positive with this 
choice). The maps were constructed by projecting the mass 
distributions within 30 h~^ kpc in the xy and xz planes. 
The colors span 4 orders of magnitude in projected density, 
with brighter colors representing higher densities. The first 
row shows the final galaxies for different feedback parame- 
ters (ec), the middle row corresponds to different energies 
per SN (-Bsn) and the lower row shows galaxies with differ- 
ent star formation efficiencies (c), as well as the no-feedback 
case. Most of our simulations show significant disks, except 
for the no-feedback experiment, E-3 (which is not very real- 
istic in terms of the energy produced per SN) and C-0.5. 

In order to quantify the properties of our disks we have 
computed the ratio between the angular momentum of each 
star in the z-direction (direction of the angular momen- 
tum of stars in the central region), j^, and the angular 
momentum expected for a circular orbit at the same ra- 
dius, jcirc: e = jz/jcirc jcirc is defined as jcirc = r ■ ^^^(rjl, 
where Vdrcir) is the circular velocity at radius r, VciTc{r) = 
yjG M{r)/r, given by the potential well of the system. A 
disk supported by rotation is characterized by a population 
of stars with e ~ 1. A component dominated by velocity dis- 
persion will show a distribution peaked near e ~ 0. In Fig. [3] 
we show histograms of e for our different simulations. The 
distribution in the no-feedback case (NF) is consistent with 
a pure spheroidal system. On the contrary, most of our sim- 
ulations with feedback produce significant disks. Different 
typical sizes and vertical extensions are found for the disks 
when the input parameters are varied (see Table[5]where the 
characteristic sizes are shown). Note that these disks coexist 
with spheroidal components which, in most cases, dominate 
the stellar mass of the system. Although the details of the 
mass distributions depend on the input parameters, it is sig- 



^ Note that our definition of e differs from that of Abadi et al. 
(2004) and Governato et al. (2007), where e is defined as the 
ratio between the angular momentum of each star in the direction 
of the total angular momentum of the galaxy and the angular 
momentum expected for a circular orbit at the same energy. 



nificant that our model produces disks for a wide range of 
parameters. 

We can use our variable e to quantify the properties 
of the disks and to separate the stars into different compo- 
nents. We have done this considering only particles in the 
main system and neglecting the satellites to avoid contam- 
ination. We decompose the stars in each simulated galaxy 
into two components which we call disk and spheroid. We 
refer to the position of the peak indicative of the disk as 
Epcak- We assume that stars with e > tpoak are part of the 
disk, and generate a symmetric distribution around epoak- 
For a given star with e > epeak, we look for a counterpart 
with e' = 2epeak — e, and with similar height above the disk 
plane, similar metallicity and similar radius. We repeat this 
for all stars with e > tpoak. This procedure avoids ad-hoc pa- 
rameters and reduces contamination by spheroid stars. As 
an example. Fig. |4] shows the edge-on and face-on maps of 
stellar surface mass density for the disk and spheroid com- 
ponents identified in E-0.7. We also show the corresponding 
velocity fields, clearly contrasting the rotation of the disk 
with the dispersion-dominated spheroid. 

We have applied our method to the runs with extended 
disk-like components. In Table [5] we show the disk and 
spheroid masses found for these cases, as well as their ratio. 
We list the values corresponding to stars within r < 2ropt 
of the central galaxy. We also list the half mass radius rdiak 
and half mass height over the disk plane /idisk for the disk 
component. From this table we can see that our experiments 
have a wide range of disk and spheroid masses. In particu- 
lar, F-0.5, F-0.9 and E-0.7 show the largest ratios between 
disk and spheroid mass. In the case of F-0.9, the disk and 
spheroidal components have the same mass, in all other cases 
the spheroid dominates in terms of mass (but may not dom- 
inate if we consider luminositieQ) . Thus, changing the input 
parameters leads to a variety of different disk sizes and scale 
lengths. 

In Table [2] we show the time when half of the mass 
of each component was assembled (rdisk and Tsphcroid). In 
the no-feedback case, the typical scale of star formation is 
very short, with Tsphcroid = 1.17 Gyr, and the spheroidal 
components in the feedback cases are also mostly old, with 
characteristic formation times below 3 Gyr. In the case of 
the disk components in the feedback experiments, their half 
mass times indicate that they are young, with typical Tdisk 
of 6 Gyr. 

Our star formation and feedback model clearly can pro- 
duce disks in a cosmological context for any choice of input 
parameters but spheroid mass dominates over disk mass 
in all cases in this particular halo, our best case being a 
1:1 disk-to-spheroid mass. In the following subsections we 
choose one of our simulated galaxies and compare its evolu- 
tion with that of the galaxy formed in the no-feedback case, 
in order to investigate in more detail how disks form and 



^ By combining the information on the age and metallicity of the 
stars, we have calculated the disk and bulge (r < 5 h~^ kpc) lu- 
minosities using the GALAXEV code (Bruzual & Chariot 2003). 
This yields a disk-to-bulge ratio of 2.4 in the R-band (and a bulge- 
to-total ratio of 0.30). This is low with respect to observational 
estimates for Milky- Way type galaxies, particularly in view of the 
fact that we neglect the effects of dust (e.g. Gadotti 2008). 
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Figure 2. Edge-on stellar surface density maps for 9 different simulations at z = 0. The upper panel shows the results for cases with 
varying feedback parameter (F-0.3, F-0.5 and F-0.9), the middle panel the series with different energy per SN {E-0.3, F-0.5 {-EsN = 1 x 10^^ 
ergs) and E-3), and the lower panel the series with different star formation efficiencies (C-0.01, F-0.5 (c = 0.1) and C-0.5), as well as the 
no- feedback simulation (NF). Note that F-0.5 belongs to all three series. 



how SN feedback affects the distributions of specific angular 
momentum, stellar age and chemical abundance. 



3.3 Evolution of the angular momentum 

Our model E-0.7 has one of the largest ratios of disk to 
spheroid mass, and has the thinnest disk. We have selected 
this simulation as an example to study in more detail the 
effects of SN feedback in galaxy evolution. The results for 
this particular simulation are qualitatively similar to those 
for all other cases where disk-like components formed. In 
order to investigate why a disk forms with feedback but not 



in the no-feedback case (NF) we compare the evolution of 
the specific angular momentum (i.e. angular momentum per 
unit mass) in the two simulations. We have computed the 
specific angular momentum for the dark matter component 
in the no-feedback and feedback simulations, using particles 
which ended up within the virial radius at z = 0. We followed 
these particles back in time, calculating their centre of mass 
and the corresponding specific angular momentum. We have 
excluded 2 = satellites which contribute an important 
fraction of the total specific angular momentum at the end, 
especially in the no-feedback case. We have also calculated 
the corresponding specific angular momentum for the cold 
gas plus stars which ended up within 2 ropt at z = 0. In this 
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Figure 3. Stellar mass fraction as a function of e = jz/jdrc for the different experiments at z = 0. The upper panel shows the results 
for the series with varying feedback parameter (F-0.3, F-0.5 and F-0.9), the middle panel the series with different energy per SN (E-0.3, 
F-0.5 (i?SN = 1 X 10^^ ergs) and E-3), and the lower panel the series with different star formation efficiencies (C-0.01, F-0.5 (c = 0.1) 
and C-0.5), as well as for the no-feedback case (NF). In the cases where a disk-like component is present, we also show the distribution 
of e for disk stars (dashed lines). In this plot, we only consider stars which belong to the main system (i.e. discounting stars in satellites). 



way, we can investigate the relation between the formation of 
the disk and the evolution of the specific angular momentum 
of the galaxy. 

In Fig. [5] we show results for the evolution of the specific 
angular momentum of the dark matter and of the cold gas 
plus stars, calculated as explained above, the left-hand panel 
shows results for the no-feedback case whereas the right- 
hand panel corresponds to E-0.7. These plots show that the 
evolution of the specific angular momentum of the dark mat- 
ter component is similar in the two cases, growing as a result 
of tidal torques at early epochs and being conserved from 
turnaround (z ~ 1.5) until z = 0. 

The specific angular momentum of the cold gas plus 
stars, however, behaves differently in the no-feedback and 
feedback experiments, in particular at late times. During 
the early stages of evolution, the specific angular momen- 



tum of this component also grows through the tidal torques 
exerted by nearby protogalaxies and it maximises at a sim- 
ilar value in the two simulations which is close to the max- 
imum value for the dark matter. However, at later times, 
the specific angular momentum of the baryonic component 
behaves differently in the two cases. In the no-feedback case 
(NF), much angular momentum is lost through dynamical 
friction, particularly through a satellite which is accreted 
onto the main halo at z ~ 1. The decrease in specific an- 
gular momentum after this time is more than an order of 
magnitude. In E-0.7, on the other hand, the cold gas and 
stars lose rather little specific angular momentum between 
z — 1 and 2 = 0. Two main factors contribute to this dif- 
ference. Firstly, in E-0.7 a significant number of young stars 
form between z = 1 and z = with high specific angular 
momentum (these stars form from high specific angular mo- 
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Figure 4. Face-on (upper panels) and edge-on (lower panels) stellar surface density maps for the final galaxy formed in E-0.7. The left- 
hand panel shows the complete stellar system, while the middle and right-hand panels show views of the disk and spheroidal components, 
respectively. These maps all have the same color-coding as Fig. [2] The arrows indicate the mean velocity field of the stars. Their lengths 
scale with velocity. 



Table 2. Mass of disk and spheroidal components (in units of 10^" h~^ Mq) and ratio between disk and spheroid masses [D / S) for the 
different simulations. We also show the disk half mass radius (r^isk) ^.nd half mass vertical scalelength (/idisk) in h~^ kpc, as well as the 
time when half of the final mass of the disk was formed (r^iak) in Gyr, as well as the corresponding time for the spheroid (Tspijcroid)- 
We also show the fraction of metals in the stellar (/™°'g) and gas (/^g') components and, for the gas component, the fraction of metals 



within 2 ropt (/^l in) ^.nd between 2 Topt and the virial radius (/" 



lit)- 



Test 



ACT 

disk 



spncroid 



D/S rdi.k 



'''disk '^spheroid /st; 



met 
ars 



irnct 
■' gas,ii 



J" met 
■/ eas.out 



NF 


- 


14.4 


- 


- 


- 


- 


1.17 


0.99 


0.01 


- 


- 


F-0.3 


1.29 


3.07 


0.42 


5.82 


1.14 


4.75 


2.95 


0.31 


0.69 


0.26 


0.74 


F-0.5 


2.46 


3.00 


0.82 


6.44 


0.84 


5.54 


2.78 


0.22 


0.78 


0.29 


0.80 


F-0.9 


3.35 


3.23 


1.04 


9.74 


0.98 


6.24 


2.33 


0.32 


0.68 


0.22 


0.78 


E-0.3 


4.85 


8.01 


0.60 


4.75 


0.66 


5.72 


2.54 


0.28 


0.72 


0.28 


0.72 


E-0.7 


3.33 


4.05 


0.82 


5.72 


0.50 


6.33 


2.49 


0.20 


0.80 


0.37 


0.63 


E-3 


- 


1.23 


- 


- 


- 


- 


2.64 


0.15 


0.85 


0.13 


0.87 


C-0.01 


2.71 


7.05 


0.39 


2.64 


0.40 


4.78 


2.83 


0.17 


0.83 


0.21 


0.79 


C-0.5 


- 


4.36 


- 


- 


- 


- 


2.38 


0.25 


0.75 


0.16 


0.84 



mentum gas which becomes cold at late times) ; and secondly, 
dynamical friction affects the system much less than in NF, 
since satellites are less massive. Our results are similar to 
those of by Zavala, Okamoto & Frenk (2007), who analysed 
two simulations with identical initial conditions but differing 
baryonic physics, one of them resulting in a disk-dominated 
galaxy, the second in a spheroidal system. 

Finally, in the right panel of Fig. [5] we show the final 



specific angular momentum of the disk and spheroidal com- 
ponents identified in this simulation at z = 0. Clearly the 
disk component has a high specific angular momentum, com- 
parable to that of the dark matter, while the spheroid formed 
in this case has a much lower specific angular momentum. 
This reflects the close relation between angular momentum 
and morphology. 
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Figure 5. Dashed lines show the specific angular momentum as a function of time for the dark matter that, at 2 = 0, lies within the 
virial radius of the system for NF (left panel) and E-0.7 (right panel). We also show with dots the specific angular momentum for the 
baryons which end up as cold gas or stars in the central 20 h^^ kpc at 2 = 0. The arrows show the specific angular momentum of disk 
and spheroid stars. 



3.4 Stellar age distributions 

Our results suggest that the formation of the disks is closely 
related to the ability to form stars at late times. In our 
no- feedback experiment, star formation occurs early, giving 
rise to a stellar spheroid sd, z — which is dominated by old 
stars. On the contrary, in E-0.7 the stellar component has 
a larger contribution from young stars. This behaviour is 
found in all simulations with feedback, as can be seen from 
Table [T] where we show the fraction of the final stellar mass 
formed between 2 = 1 and z = 0, /*^^. The no- feedback 
simulation has the lowest value for f^^^, while the largest 
values are found for the systems with the most massive disks. 

In order to analyse in more detail the relation between 
disk formation and delayed star formation, we show in Fig. [6] 
the distribution of stellar formation times for the disk and 
spheroidal components in E-0.7. The two components dif- 
fer significantly in their formation times, the spheroid being 
composed of old stars while the disk contains a younger stel- 
lar population. In particular, only 6 percent of the stars in 
the spheroidal component were formed at times later than 6 
Gyr (2 « 1). On the contrary, 53 percent of the final stellar 
mass of the disk was formed between z = 1 and 2 = 0, 24 
percent in z = [0.5,0], and 29 percent in z = [1,0.5]. This 
suggests that the ability of a system to form a disk corre- 
lates with its ability to retain enough gas to form stars at 
recent times where mergers are not as frequent as at early 
epochs. This becomes possible when SNe regulate star for- 
mation activity. We note that in our simulation without SN 
energy feedback (NF), only 10 percent of the final stellar 
mass forms after z = 1 (see Table [T|. 

The age of the disk components can also be quantified 
by the time when half of the final disk stars had already 
formed, Tdisk- The rdisk values for the different simulations 
are shown in Table[2l It can be seen that large disks typically 
have larger rdisk values. In Table[2]we also show the half mass 
formation time for the spheroidal components, Tspheroid- We 
find that, in all cases, the spheroids are mainly old. Note 
that the no-feedback case has the oldest stellar population, 
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Figure 6. Mass fraction as a function of formation time for stars 
of the disk and spheroidal components in simulation E-0.7. 



while all the feedback cases, independent of the disk they 
were able to form, show similar values of Tsphoroid- 

In order to investigate the spatial history of disk as- 
sembly, we plot in Fig. [7] the mean formation time of disk 
stars as a function of radius within the disk plane for E-0.7. 
We also show the Icr scatter and we note the times corre- 
sponding to 2: = 1 and z = 0.5. As noted above, 29 and 
24 percent of the final stellar mass of the disk was formed 
in the redshift ranges z = [1,0.5] and z = [0.5,0], respec- 
tively. Clearly the stellar disk formed from the inside-out, 
although the dispersion is high. In fact, we find that star 
formation between 2 = 1 and 2 = lies primarily in a ring- 
like structure whose radius increases with time. In contrast, 
the spheroidal component shows no trend between the age 
of stars and radius. 



3.5 Chemical abundances 

In this section we investigate the chemical properties of our 
simulated galaxies. Metals are produced in stellar interiors 
and ejected into the interstellar medium when supernova 
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Figure 7. Formation time of disk stars within 2 r^isk a-s a func- 
tion of distance within the disk plane for the feedback simulation 
E-0.7. The dotted lines give the ±cr scatter and dashed lines de- 
pict the times corresponding to ^ = 1 and z = 0.5. 



explosions take place. If feedback is neglected, most metals 
remain locked either in the stars or in the cold gas near 
star forming sites, since there is no efficient mechanism to 
transport gas and metals from the inner to the outer regions. 
Winds generated when SN feedback is included produce a 
redistribution of mass and metals, and the resulting chemical 
properties are very different from those obtained when SN 
feedback is not considered. 

Most observations of detailed chemical abundances in 
stars come from our Galaxy, and are usually quantified sep- 
arately for the different dynamical components: thin disk, 
thick disk, bulge and stellar halo. The stellar halo is found 
to be the most metal-poor component, the thin disk and 
the bulge are the most metal-rich components (e.g. Ryan & 
Nories 1991; Chiba & Beers 2000; Zoccah et al. 2003; Nord- 
strom et al. 2004). In order to study the chemical properties 
of our simulated galaxies and to see if we can reproduce these 
trends, we use the disk-spheroid decomposition explained 
above. For simplicity, we do not attempt to divide our disk 
stars into two components. In the case of the spheroidal com- 
ponent, however, we have subdivided the stars into bulge 
and stellar halo, according to distance from the centre, sep- 
arating at 5 h~^ kpc. 

In the left-hand panel of Fig. |S] we show the distribu- 
tion of iron abundance for stars in the bulge and stellar halo 
for the no-feedback case while the right-hand panel shows 
similar histograms for E-0.7 together with an abundance his- 
togram for its disk stars. In both simulations the outer part 
of the spheroid ("the stellar halo") is shifted to lower metal- 
licities than the inner part ("the bulge"). The shift is more 
pronounced in the simulation with feedback, although E-0.7 
has a lower level of enrichment overall in comparison to the 
no-feedback case. The disk and bulge components show sim- 
ilar iron abundances in E-0.7, while stars in the stellar halo 
have lower abundances, in qualitative agreement with ob- 
servation. Note that the metallicities of the inner stars (i.e. 
bulge and disk) are the result of chemical evolution in the 
main galaxy, since these components were mainly formed in- 
situ. On the contrary, stars in the stellar halo were mostly 
contributed by satellites and do not refiect the chemical evo- 
lution of the main component of the galaxy. 



The arrows in the right panel of Fig. |8] indicate how 
the distributions change for different input parameters. The 
trends are similar in all cases. An increase in tc leads to a 
shift towards higher iron abundance. Note that for higher 
ec more metals are dumped into the cold gas from which 
stars form, leading naturally to higher stellar abundances. 
On the contrary, increasing the energy per SN decreases the 
metal abundances because fewer stars are produced and the 
gas is less enriched where stars do form. As shown before, in- 
creasing the SF efficiency leads to lower stellar masses in our 
model. Consequently, metallicities are also shifted towards 
lower values in this case. 

For the gas component, the resulting metal distribution 
is affected not only by the metal production itself, but also 
by the galactic winds which transport an important frac- 
tion of the metals outwards, mixing them with the interstel- 
lar gas. For strong enough winds, metals can be lost into 
the intergalactic medium. Hence, many factors influence gas 
metallicities. In Table[2]we show the fraction of galactic met- 
als which are locked into the stellar and gas components. In 
the no-feedback case, we find that 99 percent of all metals 
are locked into stars. This result is because i) most (cold) 
gas in the central regions has been transformed into stars 
and hence most metals are locked into the stellar compo- 
nent, and ii) no winds can develop in this case so metals 
cannot be transported outwards. Consequently, gas in the 
outer regions (where star formation is less efficient) is not 
enriched. On the contrary, when our feedback model is in- 
cluded, the fraction of metals locked into stars significantly 
decreases, being typically of the order of 20 — 30 percent. 
Because of the efficiency of SN feedback in blowing galactic 
winds and regulating star formation, most of the metals are 
now in the gas component. 

In order to quantify the importance of winds in dis- 
tributing metals through the gas component, the last two 
columns of Table [2] show the fraction of metals in the gas 
located in the inner regions (r < 2 ropt, fglT.in) ^i^d in the 
outskirts of the systems (2 ropt < r < r2oo, fglT.out)- In the 
no-feedback case, only 1 percent of the metals are in the 
gas. For our feedback cases we find typically that 70 percent 
of the gaseous metals are outside the innermost region. In 
Fig. |9]we show the oxygen profiles for the gas components in 
simulations NF and E-0.7. SN feedback strongly affects the 
chemical distributions. If no-feedback is included, the gas is 
enriched only in the very central regions. Including SN feed- 
back triggers a redistribution of mass and metals through 
galactic winds and fountains, giving the gas component a 
much higher level of enrichment out to large radii. A lin- 
ear fit to this metallicity profile gives a slope of —0.048 dex 
kpc~^ and a zero-point of 8.77 dex, consistent with the ob- 
served values in real disk galaxies (e.g. Zaritsky et al. 1994). 

We have shown that our star formation and feedback 
model produces galaxies with well differentiated stellar dy- 
namical components which also differ in their chemical abun- 
dances. In addition, it can produce the correct gas metallic- 
ity profiles, although our chemical model is still very simple. 
Galactic winds produced by SN feedback redistribute metals 
within galactic haloes. This is an important success of the 
model. It opens up the possibility of studying the detailed 
distribution of metals within galaxies forming in their full 
cosmological context, as well as exploring the enrichment of 
the intergalactic medium. 
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Figure 8. Mass fraction as a function of iron abundance for stars in our no-feedback simulation (NF, left-hand panel), and in E-0.7 
(right-hand panel panel). We show separately the results for stars in the disk (for E-0.7), the bulge and the stellar halo. Vertical lines 
indicate [Fe/H]=-1 and [Fe/H]=0. 
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Figure 10. SFRs for simulations of a ~ 10 h Mq virial mass system run without (DNF) and with (DE-0.7) energy feedback, as well 
as for TE-0.7, a ~ 10^" h~^ Mq simulation including feedback. To facilitate comparison, the SFRs arc normalized to the scale factor F. 
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Figure 9. Oxygen abundance for the gas component as a func- 
tion of radius projected onto the disk plane for our no- feedback 
simulation (NF) and for the feedback case E-0.7. The error bars 
correspond to the standard deviation around the mean. 



4 DEPENDENCE ON GALAXY MASS 

In this Section we investigate how the effects of SN feed- 
back change with galaxy mass. In Scannapieco et al. (2006) 
we showed that the star formation is affected differently by 
SN feedback in systems of differing mass. We found that for 
large systems (virial masses of order Mvir ~ 10^^ /i^^Mq), 
SN feedback reduces and prolongs star formation with re- 
spect to simulations where it is not included. In smaller sys- 
tems (Mvir < 10^° /i~^Mq), star formation not only de- 
creases but its character also changes, becoming much more 
bursty. Small systems are very strongly affected by SN feed- 
back, developing galactic winds that are able to expell most 
of their baryonic mass. In this Section, we again test the 
effects of SN feedback on smaller mass systems in order to 
see if the results found in Scannapieco et al. (2006) are valid 
for galaxies forming in their proper cosmological context. 

We have scaled down the initial conditions for simula- 
tion E-0.7 by factors of F = 10~^ and 10^'^ in mass and 
r^' ^ in linear scale to generate two smaller systems of final 
virial masses ~ 10^ h~^ Mq (DE-0.7) and ~ 10^° h'^ Mq 
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(TE-0.7). These simulations were run with the same input 
physics as E-0.7. These initial conditions allow us to compare 
simulations of galaxies of very different mass but identical 
assembly history. Note however they may have biased envi- 
ronments and assembly histories compared to real dwarfs. 
We have also run a simulation of the ~ 10 h~^ M© virial 
mass galaxy without energy feedback from SN (DNF). This 
is a low mass counterpart of NF. In Table [3] we indicate the 
scale factor F, as well as the dark matter, gas and stellar 
masses of the final z — galaxies renormalised by the scale 
factor for these simulations. This is equivalent to expressing 
masses in units of lO'' /I'^M© for DNF and DE-0.7, and in 
units of 10* /i"^Mq for TE-0.7. We also list the gas and 
baryon fractions for the final galaxies in these simulations. 
All these quantities are computed within the virial radius. 
For comparison, in the lower row we repeat the results for 
E-0.7 from Table [T] In this way, the results shown in this 
table can be directly compared; any difference reflects the 
way SN feedback varies with virial mass. 

In the left-hand panel of Fig. [10] we show SFRs renor- 
malized by the scale factor for the ~ 10^ h~^ M© virial 
mass simulations with and without SN feedback. SN feed- 
back clearly has a dramatic effect on the SFR of such small 
galaxies. When SN feedback is not included (DNF), the SFR 
behaves similarly to that in NF (a 10^ times more massive 
galaxy), leading to a final stellar mass of 1.88 x 10* h~^ Mq. 
Within the virial radius, the final mass fraction in stars is 
0.09 which is comparable to that for NF (0.11). The SFR in 
DE-0.7 behaves quite differently. It is very low at all times; 
SN feedback produces a decrease in SFR of about 2 orders 
of magnitude compared to the no-feedback case. The final 
stellar mass formed is only 10*' h~^ M© (~ 1200 star parti- 
cles) in this case. The amount of gas within the virial radius 
is also smaU, Mgas = 3.7 x lO'' H'^Mq. This is the result 
of the violent winds that develop in this galaxy due to the 
shallow potential well. As a consequence, most gas fails to 
condense and form stars. From the last two columns of Ta- 
ble |31 we can see that the baryons in DE-0.7 represent only 
2 percent of the final mass in the galaxy, as opposed to the 
case without SN feedback where the final baryon fraction is 
10 percent. 

In the right-hand panel of Fig. 1101 we show SFRs renor- 
malised by the scale factor for DE-0.7, TE-0.7 and E-0.7 
which differ only in virial mass and were all run with the 
same SN feedback model with identical input parameters. 
The effects of SN feedback increase very substantially with 
decreasing virial mass. Table [3] shows that more massive 
galaxies at z = tend to have larger gas and stellar mass 
fractions (AZ/F). This is because more of their gas is able to 
cool and condense to build up the stellar component, and 
less gas is lost in galactic winds. For these reasons, the larger 
the galaxy, the lower the final gas fraction (gas over bary- 
onic mass), and the larger the baryon fraction (baryonic over 
total mass). 

We note that poor resolution may affect our results for 
the smallest galaxy. In Scannapieco et al. (2006) we carried 
out resolution tests of our feedback implementation. This 
analysis suggested that small galaxies must be represented 
by many particles (at least 40000) to get convergent results 
because their final state is heavily infiuenced by details of 
the interplay between heating and cooling. Here our galax- 
ies form through both merging and accretion and it is not 



Table 3. List of properties of the simulations with final virial 
mass ~ 10^ h-^ M© and ~ 10^" h'^ Mq. We show the scale 
factor r relative to E-0.7, the dark matter, gas and stellar masses 
renormalised by the scale factor F in units of 10^" h~^M.Q. We 
also show the gas and baryonic fractions within the virial radius. 
For comparison, we show the values for E-0.7 from Table [T] 
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clear how much poor resolution at high redshift affects our 
results at 2 = 0. In the simulations presented in this paper, 
all the large galaxies, as well the smaller object DNF, are 
represented at 2 = by more than 100000 dark matter and 
70000 baryonic particles. Only in runs DE-0.7 and TE-0.7 is 
the number of baryonic particles ( 17000 and 40000, respec- 
tively) so small that better resolution is probably needed to 
obtain full convergence. However, although the final proper- 
ties of the systems may slightly differ at higher resolution, 
our general conclusions about the dependence of SN feed- 
back effects on virial mass should still be valid. 

The results obtained in this section prove that our 
model is able to reproduce the expected dependence of SN 
feedback on virial mass without changing the relevant phys- 
ical parameters. Thus, our model is well suited for studying 
the cosmological growth of structure where large systems are 
assembled through mergers of smaller structures and sys- 
tems form simultaneously over a wide range of scales. 



5 CONCLUSIONS 

We have used a series of simulations of the formation of 
a galaxy of virial mass ~ 10^^ /j-i M© from cosmologi- 
cal initial conditions to study how SN feedback affects the 
formation of disks. Our simulations have sufficient resolu- 
tion to resolve the internal structure of the galaxy as well 
as large-scale processes such as mergers, interactions, infall, 
tidal stripping, etc. Our study has focussed on the perfor- 
mance of our SN feedback scheme for realistic initial con- 
ditions and on the dependence of the results on the input 
parameters of the feedback model, as well as on the mass of 
the galaxy. In the following we summarize our main results. 

1) Our model for SN feedack is able to regulate star 
formation activity by heating and pressurizing gas near star- 
forming sites. This is refiected in a bursty behaviour for 
the SFR when feedback is included, particularly in low-mass 
systems. Different choices for feedback efficiencies produce 
variations in the star formation levels and in the detailed 
shape of the SFR, although qualitatively all our simulations 
show similar behaviour. When compared with a simulation 
where SN feedback is neglected, the reduction in the final 
stellar mass formed can be as high as 70 percent even for a 
Milky Way mass object. 

2) Efficient regulation of star formation leads to an in- 
crease in the final gas fraction, when compared to a simula- 
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tion where SN feedback is not included. This, together with 
substantial mass-loaded galactic winds powered by the SN 
input, leads to galaxies that, at 2 = 0, have lower baryon 
fractions. 

3) Galaxy morphology is strongly affected by SN feed- 
back. Large disks can be formed and can survive until z = Q 
if energy injection is allowed. Our results indicate that disks 
can, however, be formed for a wide range of input parame- 
ters. The stellar mass in our galaxies is typically dominated 
by spheroidal components. In the best case, we found equal 
masses in the disk and the spheroid. Note, however, that the 
disk is significantly more luminous since it is significantly 
younger. We find the final characteristic size and thickness 
of disks to depend significantly on the input parameters. 
Adopting too extreme and unrealistic assumptions for the 
energy per SN or for the efficiency results in failure to form 
a disk. Similarly, a simulation with no-feedback also leads 
to a pure spheroidal stellar component. 

4) Disks are rotationally supported, have high specific 
angular momentum and have a substantial number of young 
stars. Our disks have typical half-mass formation times of 
~ 6 Gyr. In contrast, spheroids are old, with half-mass for- 
mation times < 3 Gyr. In a simulation without feedback 
much of the specific angular momentum of the cold compo- 
nent (the cold gas and stars) is lost to dynamical friction 
during the assembly. As a result, no disk forms. 

5) Our disks and spheroids have different chemical 
properties as a result of their different formation histo- 
ries. In particular, the stellar halo in our simulations is the 
most metal-poor component, in agreement with observation, 
whereas the disk and bulge components have similar near- 
solar metallicities. Without SN feedback this trend is not 
recovered and all components show similar, overly high iron 
abundances. Gas metallicities are also affected by changes 
in star formation histories and by galactic winds when SN 
feedback is included. Gas metallicity profiles consistent with 
observations are only found when feedback is included. 

6) In our model, the efi'ects of SN feedback strongly de- 
pend on virial mass. We have shown this by scaling down one 
of our simulations by factors of 10^ and 10^ in mass. Feed- 
back has a much stronger efi'ect in smaller systems, where 
violent galactic winds eject most of the baryonic material 
and prevent gas from collapsing to make stars. SFRs are 
low and bursty in such systems. Such small galaxies have 
lower baryon fractions at 2; = than larger systems. 

Our findings show that our model is able to produce 
realistic disk systems from cosmological initial conditions. 
This results from the self-regulation of the star formation 
process and the heating due to feedback which prevents ex- 
cessive early condensation of gas into the central regions of 
dark haloes and the consequent later loss of angular mo- 
mentum. An important success is the fact that disks are 
obtained for a variety of input parameters, with realistic 
radial and vertical sizes, and metallicity structures. In all 
cases we find that the spheroids dominate the stellar mass, 
suggesting that further refinements are needed if we are to 
understand the formation of late-type spirals. 
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